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In addition to the well known ferromagnetism, double perovskites are also expected to exhibit 
antiferromagnetic (AF) order driven by electron delocalisation. This has been seen in model Hamil- 
tonian studies and confirmed via ab initio calculations. The AF phases should occur, for example, 
on sufficient electron doping of materials like Sr2FeMo06 (SFMO) via La substitution for Sr. Clear 
experimental indication of such AF order is limited, possibly because of increase in antisite disorder 
with La doping on SFMO, although intriguing signatures of non ferromagnetic behaviour are seen. 
We study the survival of electronically driven antiferromagnetism in the presence of spatially cor- 
related antisite disorder and extract the signals in magnetism and transport. We discover that A 
and G type AF order, that is predicted in the clean limit, is actually suppressed less strongly than 
ferromagnetism by antisite disorder. The AF phases are metallic, and, remarkably, more conducting 
than the ferromagnet for similar antisite disorder. We also highlight the phase coexistence window 
that connects the ferromagnetic regime to the A type antiferromagnetic phase. 



I. INTRODUCTION 

The double perovskites 1 3 (DP) are materials of the 
form ABOsAB'03. The B and B' are usually transition 
metal ions, while A is a rare earth or alkaline earth. The 
itinerant electrons in this system are coupled strongly to 
the transition metal magnetic moments, which act as core 
spins, and the magnetic order is driven 4 5 by minimisa- 
tion of the electronic energy, rather than a short range 
interaction between the moments. Much of the excite- 
ment in the double perovskites has been due to the high 
ferromagnetic T c , noticeable magnetoresistance, and the 
possibility of spin polarised conduction in materials like 
Sr 2 FeMo0 6 (SFMO). 

While the ferromagnetism is certainly useful, one may 
wonder about the occurence of other kinds of magnetic 
order in these materials. For example, even 'simple per- 
ovskite' transition metal oxides, the cuprates, mangan- 
ites, or cobaltates, have a rich phase diagram 6 , with a 
strong dependence on the doping level. The manganites, 
for instance, exhibit not just ferromagnetism, but also 
'CE' magnetic order and A, C, and G type antiferromag- 
netic (AF) phases-^, depending on the hole doping level. 
It is interesting to explore if non-ferromagnetic ordered 
states are possible in the double perovskites as well. 

A study of the model Hamiltonian of these materials 
shows that at low electron density a ferromagnetic (FM) 
alignment of the core spins is favoured since it leads to 
the maximum bandwith. However, at sufficiently large 
band filling, antiferromagnetic (AF) states with A or 
G type order (in two dimensions) successively become 
favourecPE23 While these spin configurations lead to 
smaller electronic bandwidth they have a higher density 
of band edge states compared to the ferromagnet. 

We have studied these magnetic states in two 
dimensions 8 and three dimensions 11 elsewhere. Ab ini- 
tio calculations using simple collinear arrangement of the 
'core spins' predict similar results^^ in the context of 
real materials. Nevertheless, electronically driven AF or- 



der continues to be elusive experimentally. One can an- 
ticipate two reasons: (i) the AF phase occurs at (high) 
electron densities which have not been probed yet, or 
(ii) increasing electron density leads to a rapid growth in 
antisite disorder (ASD), strongly suppressing any signa- 
ture of long range AF order. 

The theoretical effort till now has focused on AF order 
in the ideal structurally ordered background, where the 
B and B' ions of A2BB'06 alternate along each axis. The 
situation in the real material is far from ideal, and the 
B, B' alternation is interrupted by an antiphase bound- 
ary (APB/^E^ involving BB or B'B' nearest neighbours. 
This leads to a pattern of structural domains. The ASD 
not only destroys periodicity but also brings into play 
an additional AF superexchange coupling when two B 
ions (like Fe) adjoin each other. If this superexchange 
scale is sufficiently large, the BB antiphase boundary 
also acts as a magnetic domain walP^ (MDW). The ef- 
fect of such AF superexchange in the double perovskite 
ferromagnets is well known. There it leads to magnetic 
domain formation, suppressing the magnetisation, and 
large enhancement of the resistivit y 14 * 15 *. The interplay 
of antisite disorder and AF superexchange with itinerant 
electron antiferromagnetism is unexplored. 

This paper provides a systematic exploration of the 
effect of increasing antisite disorder on the AF phases 
in a two dimensional (2D) double perovskite model. We 
work in 2D for ease of visualisation and to access large 
system size. 

The phase diagram mapping out the occurence of AF 
phases in the 'clean' limit has been established earliei^. 
Here we focus on a couple of electron densities, one each 
in the 'A type' and 'G type' window, respectively, and 
study the magnetic order and transport for increasing 
antisite disorder. Our principal results are the follow- 
ing, (i) The suppression of magnetic order with ASD is 
slower in the AF phases compared to the ferromagnet, 
and even with ~ 30% mislocation there are clear signa- 
tures of long range order, (ii) The T c is not significantly 
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affected, until strong disorder, although the 'transition' is 
broadened, (iii) Both AF phases are metallic, and the in- 
crease in their resistivity with ASD is weaker than in the 
ferromagnet. (iv) The phase separation regime between 
the FM and A type phase is converted to a phase coexis- 
tence window in the presence of moderate ASD, and one 
observes the occurence, locally, of both kinds of order in 
the same sample. 

The paper is organised as follows. We first quickly re- 
view the results on AF phases obtained in model studies 
and ah initio calculations, and follow it up with a survey 
of the few experiments in this regime. We then discuss 
the model and method. This is followed by a discussion of 
our results on the AF order in the background of antisite 
disorder, with focus on the spatial correlations and the 
transport properties of the system. We conclude with 
a discussion of the phase coexistence region that inter- 
venes between the ferromagnetic and A type phase, and 
is likely to be of relevance in La doped SFMO. 



II. STUDIES ON AF ORDER 



Early studies using model Hamiltonians for double per- 
ovskites had observed the instability 5 of the ferromag- 
netic state, without exploring the competing phase that 
emerges. A subsequent variational study 16 did iden- 
tify non-ferromagnetic phases. More recent studies using 
both simple models^ and realistic DFT calculations^^ 
indicate that the ferromagnet becomes unstable to 
an A type phase on increasing electron density. In 
S^-aXa^FeMoOe, for example, this is expected to hap- 
pen for x > 1. The DFT studies have employed su- 
percells for a few commensurate doping levels. Using a 
three band model Hamiltonian, with parameters inferred 
from the DFT, the authors have explorecP a more con- 
tinuous variation of La doping level and confirmed the 
DFT trends. The crossover to a non FM ordered state 
is, therefore, not an artifact of a single band model or 
two dimensionality that earlier studies employed. As for 
the effect of ASD on non ferromagnetic phases, we are 
aware of only one study involving uncorrelated antisite 
defects 16 . It is more focused on the doping dependence, 
and explores mainly the magnetism, but the trends are 
consistent with what we observe here. 

Samples have indee d been synthesised with large La 
doping on SFMO 1718 . The main observations are (i) a 
suppression^ of the low field magnetisation with increas- 
ing x, and (ii) a large difference^ between the field cooled 
(FC) and zero field cooled (ZFC) response. There is 
unfortunately no detailed understanding of the ASD in 
these samples yet (for example via XAFS), or data on 
resistivity and magnetoresistance. 



III. MODEL AND METHOD 

The model we study has been presented earlier in the 
context of the ferromagnetic phase, so we only provide a 
brief description for completeness. 

The information about the B, B' structural motif is 
encoded in a binary variable 77^ with rji = 1 for B sites 
and r]i = for B' sites. In the structurally ordered DP 
the rji alternate along each axis. We will consider pro- 
gressively 'disordered' configurations, generated through 
an annealing process. For any specified {77} background 
the model has the form: 

H = e B ^2vifLfio- + e B'^2(l ~ r)i)m\ G m i(J 

+JAF^2viVjSi ■ Sj (1) 

(id) 

f is the electron operator on the magnetic B site 
and m is the operator on the non-magnetic B' site. 
€b and cb' are onsite energies, at the B and B' sites 
respectively. e# — e#/ is a 'charge transfer' energy. 

is the electron hopping term: —t\ ? v ViVjfLfjo- 
-*2 £<i,j>(l - Vi)0- - Vj)mtm ja -t 3 T,(i,j)*(Vi + Vj ~ 
2r]ir]j)(fJ a mj a + h.c). The t's are all nearest neighbour 
hopping amplitudes, and for simplicity we set £]_ = t 2 = 
ts = t here. is the core spin on the site R^, with 
I Si I = 1. J is the Hund's coupling on the B sites, and we 
use J/t ^> 1. When the up-spin core levels are fully filled, 
e.g., Fe anc ^ e #T m SFMO, the conduction electron is 
forced to be antiparallel to the core spin. We have used 
J > to model this situation. For the present study we 
have set 19 the effective level difference e# — J/2 — = 0. 
The chemical potential \i is used to control the electron 
density, and TV is the total electron density operator. We 
study one point in the A type window, and another in 
the G type window. We also explore the transition from 
FM to A type order. 

When two magnetic atoms can be on neighbouring 
sites, we also have to consider the antiferromagnetic cou- 
pling Jaf between nearest neighbour B sites. 

It has been well established now that the antisite dis- 
order does not involve 'random' exchange of B and B' 
ions, but follows a correlated pattemP2E] j± periodic 
B-B' pattern is interrupted by a line-like defect, where 
BB or B'B' adjoin each other, and beyond this boundary 
one obtains another domain but with a phase slippage 
with respect to the first one. The line defect is an an- 
tiphase boundary. If a fraction x of atoms are 'mislo- 
cated' with respect to the ideal ordered structure, these 
atoms themselves are organised into domains, and do not 
act as simple 'point defects'. 

We generate the correlated patterns by studying a lat- 
tice gas model 20 , with poor annealing to prevent long 
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average, (1/Acj) a(uj)duj of the optical conductivity. 
This is averaged over thermal configurations and disor- 
der, as appropriate. We use Aoo = O.Obt. 






IV. RESULTS 




FIG. 1: Top: The three magnetic phases in the structurally 
ordered 2D double perovskite model. Left: ferromagnetic 
(FM), center: A type antiferromagnet, right: G type anti- 
ferromagnet. These occur with increasing electron density. 
The moments are on the B sites. We have not shown the in- 
duced moments on the B' sites. Bottom: Electronic density 
of states for the ferromagnetic, A and G type ordered phases 
in the structurally ordered background. 



range order. Briefly, the {rji} configuration arise from: 



H eff { V } = -Vj2vi(l 



■Vj) 



(2) 



The r] are coupled only between nearest neighbour sites. 
The ground state in this model corresponds to 77 = 
1,0,1,0,.. along each axis, i.e, B, B', B, B'.. V > 
is a measure of the ordering tendency, and sets the tem- 
perature for long range B-B' order if the system were 
allowed to equilibriate. 

We used different 'annealing time', r ann , and annealing 
temperature, T ann , to generate the ASD configurations. 
The ASD structures are characterised in terms of the fol- 
lowing indicators: (i) The structural 'order parameter' 
S = 1 — 2x, where x is the fraction of B (or B') atoms 
that are on the wrong sublattice. (ii) The degree of short 
range order, characterised by the probability, p, of having 
nearest neighbour pairs that are B-B'. (Hi) The corre- 
lation length £ in these structures, computed from the 
width of the ordering peak. 

For a given {rji} configuration, generated at some 
T ann and r ann , we solve for the magnetic properties, 
and electronic properties averaged over equilibrium mag- 
netic configurations, via a cluster based Monte Carlo 
technique 23 . Electronic properties are calculated after 
equilibration by diagonalizing the full system. 

The electronic conductivity is calculated using the 
Kubo formula, computing the matrix elements of the cur- 
rent operator. The 'dc conductivity' is the low frequency 



A. The clean limit 

In the absence of antisite disorder there are no AF su- 
perexchange interactions in the system, and the magnetic 
order is decided by minimisation of the electronic energy. 
We reproduce the model for the ordered DP, below, for 
ease of reference: 



H rd = 6 B ^2 fief™ + € B' ^2 m L m iv ^ 
ieB,a ieB'a 



<ij>cr 



ieA,a{3 



For J/t ^> 1 this model supports three collinear phases 
(in 2D). The FM state gives way to A type (line like) 
order with increasing electron density, and finally to a 
G type state. These have been discussed earlier^, we 
reproduce the magnetic configurations, top row in Fig.l, 
and the electronic density of states below. We have set 
(e B - J/2) - e B > = 0. 

The FM state is preferred at low n, since it has the 
largest bandwidth. The A type state has lower band- 
width, but with large density of states near the band 
edge. The FM becomes unstable at n ^ 0.45. The A 
type state is stable for n > 0.58, and between these we 
have a phase separation (PS) window. Similarly the A 
to G transition involves a PS window. The PS windows 
narrow with increasing T and vanish as T —> T c . The 
thermal transitions and the PS windows are shown in 
Fig.2. 

Broadly, the task would be to extend this phase dia- 
gram to finite antisite disorder. Instead of attempting 
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FIG. 2: Colour online: Phase diagram for the non disordered 
double perovskite. We only show the region n — [0, 1]. From 
n — 1 to n — 2 one populates the non-dispersive B r level, 
and the magnetic state is G type. The n — [2,3] window 
is a symmetric version of the n — [1,2] region. The regions 
between the phases indicate phase separation. The results are 
obtained via MC on a 40 x 40 lattice. 
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FIG. 3: Colour online: Antiphase domains and correspond- 
ing antiferromagnetic phases. The left column shows the do- 
main pattern in the ASD background, with structural order 
parameter S = 0.98, 0.76, 0.50, 0.08 as we move from top 
to bottom. The middle column shows the A type AF phase 
on this structural motif, while the right column shows the G 
type phase. The magnetic correlations are characterised via 
the overlap factor gi = So. Si, where So is the left lower corner 
spin in the lattice. 



to map out the disorder dependence at all densities we 
choose two representative densities, n ~ 0.65 in the A 
type window, and n ~ 0.95 in the G type region, to clar- 
ify the impact of disorder. We also explore the effect of 
ASD on the phase separation window since it would be 
encountered in any attempt to electron dope the ferro- 
magnet. 



B. Disorder configurations 

We study four families, with progressively increasing 
antisite disorder. Each family arises from annealing the 
lattice gas model for some time r ann at temperature T ann . 
A representative configuration from each family is shown 
in the first column in Fig. 3. They have a fraction of 
mislocated sites: x = 0.01, 0.12, 0.25, 0.46. The frac- 
tion of course varies somewhat from sample to sample 
within each family. Our 'disorder average' for magnetic 
and electronic properties is performed typically over 10 
configurations within each family. 

The structural 'order parameter' S = 1 — 2x is 



0.98, 0.76, 0.50, 0.08 for these four families. Since the 
configurations emerge from an annealing process, the 
spatial correlations in the structure are much stronger 
than the S value would suggest. In the absence of 
information about spatial correlations, the probability 
p of having a BB' nearest neighbour bond is p = 
x 2 + (1 - x) 2 = (1/2)(1 + S 2 ). The two terms, x 2 
and (1 — x) 2 , arise from having both atoms on the 
'wrong' sublattice or 'right' sublattice respectively. For 
the S values in our configurations these numbers are ~ 
0.98, 0.79, 0.63, 0.50. In our case the structures are spa- 
tially correlated, so the probability of BB' nearest neigh- 
bours is much higher than the estimates, above, from an 
'uncorrected distribution'. By analysing our configura- 
tions we obtain: p ~ 0.98, 0.97, 0.95, 0.86, indicating a 
high degree of short range order. We also estimated the 
'correlation length' £ associated with the domain struc- 
tures, from a Lorenztian fit to the B-B' structure factor 
of the form Sbb 1 (q) ~ 7((<fe - ^? + (% ~ ^? + T 2 )- 
This yields £ - 6.6, 5.9, 4.8, 3.6. 

Overall, the disorder in these systems should be char- 
acterised in terms of two variables: S, which is a gross 
measure of order, and p (or £) which quantify short range 
correlation. In general, physical properties would depend 
on both of these and not simply S. 



C. Magnetic order with antisite disorder 

In the FM case it is simple to see that the presence of 
AF superexchange at the APB would tend to align spins 
in opposite directions across an antiphase boundary. The 
system breaks up into up and down spin domains. Sup- 
pose the 'up spin' domains correspond to the correctly 
located sites and are the 'majority'. The net magnetisa- 
tion is proportional to the volume difference between the 
correctly located and 'mislocated' regions. If the degree 
of mislocation is x, then the normalised magnetisation 
M = (1 — x) — x = 1 — 2x = S. In the magnetic struc- 
ture factor £>(q), the FM peak is at Q F = {0,0}. By 
definition D(Qp) = M 2 , and from the domain argument 
this is simply S' 2 . This dependence is well established 
experimentall^EH, and also observed by us^l 

For the A type phase, studied on the same antisite 
structures as the FM, the nature of local magnetic or- 
der is more subtle. For the clean A-AF the order is 
at two possible pairs of wavevectors, Q^i = {7r/2,7r/2} 
and Qa2 = {37r/2, 37r/2}, or Q^3 = {7r/2,37r/2} and 
Qa4 = {37r/2, 7r/2}. We have set the lattice spacing 
ao = 1 on the DP lattice. The two sets arise due to the 
two possible diagonals along which the FM stripes can or- 
der. Within each set there are two Q values because half 
the sites in the DP lattice are non-magnetic, and the spin 
field has to have nodes there. In the clean system either 
{Qai?Qa2} or {Qa3,Qa4} are picked. In a disordered 
system all four can show up, as in the middle column 
bottom row in Fig. 3. For the G type phase the order is 
at the single pair: Qgi = {7r, 0} and Qg2 = {0,7r}. 
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FIG. 4: Colour online: Magnetic order in the A type and 
G type phases with ASD. The results are based on MC on 
40 x 40 systems, and averaged typically over 10 configurations 
for each value of S. 

The middle and right columns in Fig. 3 show A type 
and G type order, respectively, for progressively increas- 
ing ASD. Fig. 4 quantifies the suppression of the ordering 
peak in -D(q), after disorder averaging over copies with 
roughly fixed degree of mislocation. 

We will analyse the structure factor, D(q), in terms 
of the domain pattern. D(q) is related to the Fourier 
transform of the spin configuration: 

0(q) = ^l/(q)l 2 

/(q) = £s i e i * R < 

Ri 

where V is the total volume of the system. 

For collinear order, where the spin projection is only 
on the z axis, the spin vector can be replaced by Si. If 
we imagine the spin configuration to be broken up into 
domains, indexed by a label a, say, then, for collinear 
phases, /(q) = zf(q) and: 

/(q) = EE^ Jlf = E^( , J) ( 3 ) 

a Hf a 

where the sum runs over the domains, and the rf are 
coordinates within a domain. This shows that /(q) gets 
additive contributions from various domains, with phase 
factors that we will soon clarify. The formulation above 
holds as long as (a) there is no significant non collinearity 
and (b) £ ^> 1, i.e., we can ignore 'interfacial' spins which 
may be hard to assign to any particular domain. 

1. A type order 

For weak to moderate ASD we observe that the sys- 
tem prefers FM stripes along any single diagonal, albeit 
with phase slippage between the stripes to accommodate 
the effect of Jaf- This is true for S = 0.98, 0.76 and 
0.50, the top three rows in the middle column in Fig. 3. 



Domains which are translated with respect to the ref- 
erence domain have relative displacement 5r a = xao or 
yao. The order within all domains is similar. So, the con- 
tribution of each domain at the ordering wavevector, Q, 
will be proportional to the domain volume, and involve 
a phase factor: 

/a(Q) = V r a e < ^ r -/o(Q) 

where /o(Q) is the normalised reflection in the perfectly 
ordered system. For the Q = Q^i = {7r/2,7r/2} peak, 
the phase factor is e t7V ^ 2 irrespective of whether the do- 
main is x displaced or y displaced. So, all the 'mislocated' 
sites, grouped into domains, contribute V m i S e t7r / 2 fo(Q), 
where Vmis is the total volume of mislocated regions. 

Adding the contribution from the majority domains, 
which are undisplaced, we obtain: 

/(Qai) = ((V ~ Vmis) + V mzs e'^ 2 )f (Q A1 ) 

Remembering that V m i S /V = x, the volume normalised 
structure factor peak is 

D(Qai) = 1(1 -x) + e"' 2 x\ 2 = 1(1 + S 2 ) 

This is roughly consistent with the S dependence of the 
T — )> structure factor in Fig. 4. (a). It is distinctly slower 
than the suppression of order in the ferromagnet, where 
D(Q F ) ~ S 2 . 

At larger ASD however, the system has short stripes 
oriented along both diagonals, see bottom row middle 
panel in Fig. 3. These domains have magnetic peaks at 
Qa3, Qa4 and not at Qai, Qa2 Even assuming that the 
'majority' domains all contribute (1 — #)/o(Qai)> we no- 
tice that the mislocated regions require classification into 
two groups: those contributing at Qai, Qa2 with volume 
fraction y, say, and those at Qa3, Qa4 with volume frac- 
tion x — y. In that case the peak at Q^i would be 

D(Q A1 ) = \(l-x) + e^ 2 y\ 2 

The [x — y) fraction makes no contribution to the peak, 
and that weight is 'lost'. Notice that x > y > 0, and it is 
not possible to write D(Qai) purely in terms of S. We 
could write the expressions for the structure factor at the 
other three Q as well, and they will all depend on both x 
and y. This is a general feature of magnetic states where 
the order can locally pick out different orientations. 

We can make some headway in the strong disorder 
limit, x = 1/2, S = 0, by assuming that there are four 
kinds of domains, with roughly equal area. There would 
be two families of {Qai?Qa2} domains, each with 1/4 
the system volume, and a relative phase shift tt/2. Simi- 
larly there would be two families of {Qa3, Qaa} domains, 
each with volume 1/4, and relative phase shift it/ 2. In 
this case the Q^i peak, for example, in D(Q) will be: 

D(Q A1 ) = |l/4 + e-/ 2 /4| 2 = l/8 

This is not very far from D(Qai) ~ 0.1 that we obtain 
from the Monte Carlo. Other peaks, at Qa2,Qa3> 
would have similar magnitude. 
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2. G type order 

G type order occurs at the combination {Qgi> Qg2} • 
{{0,7r}, {7r,0}}. As before the relative displacement of 
the domains can be only xao or yao. Suppose we are 
computing the structure factor at Qgi then all domains 
will contribute, but with following phase factors: zero if 
the domain is not mislocated (5r = 0), zero again if the 
domain is x displaced, and e l7T = —1 if the domain is y 
displaced. 

In two domain systems, as in the second row in Fig. 3, 
the mislocated domain is either x displaced or y dis- 
placed. For copies with x displacement the contribution 
at Qgi will be |(1 — x) + x\ 2 = 1, while for y displace- 
ment it will be |(1 - x) + e i7r x\ 2 = (1 - 2x) 2 = S 2 . Av- 
eraging the structure factor over copies would lead to 
D(Q G i) = (1/2)(1 + S 2 ). This is roughly what we ob- 
serve in our Fig.4.(b) at T = 0. 

In large systems, where there will be many domains, 
we can assume that half the mislocated domains are x 
displaced and half y displaced. In that case the structure 
factor would be 

D(Q G i) = 1(1 -x) + x/2 + e^x/2\ 2 

Using 1 — 2x = S this leads to 

D(Q G1 ) = (1/4)(1 + S) 2 

For S — >• this gives 0.25, not far from ~ 0.20 that we 
obtain from our configurations. 

This reveals that for both A and G type order even 
when half the sites are mislocated, i. e, one has maximal 
antisite disorder, there is a surviving peak in the struc- 
ture factor. All these of course assumed that the struc- 
tural pattern had a high degree of spatial correlation so 
that one can meaningfully talk of domains. We should 
have 1 — p <C 1, or structural correlation length £ ^> a$. 
If the structures were fragmented to a random alloy then 
the results above would not hold. We have checked this 
explicitly. 
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FIG. 5: Colour online: Resistivity in the A type and G type 
phases with ASD. The results for each S are averaged over 10 
realisations of disorder. 



FIG. 6: Colour online: Comparison of the resistivity in the 
FM, A and G type correlated phases in samples with the 
lowest degree of order S = 0.08. 



D. Transport with antisite disorder 

All the three phases, FM, A type, and G type, in the 
2D double perovskite model are metallic in the clean 
limit. The electronic states are extended, and there is 
a finite density of states at the Fermi level. In the ab- 
sence of antisite disorder the resistivity, p(T), in all three 
have similar temperature dependence. The resistivity in- 
creases rapidly as T increases towards T c , Fig. 5 and our 
earlier worlJ 2 ^] on the FM, and 'saturates' at high tem- 
perature. A Fisher-Langer type 2 -^ phenomenology can 
qualitatively describe the transport. 

Weak disorder leads to an increase in the residual re- 
sistivity of the A and G type phases, see Fig. 5, as ob- 
served earlier for the FM. The sharp resistive transition 
observed in the clean limit (S = 0.98) is also gradually 
broadened with increasing disorder. There is however a 
key difference with respect to the FM phase when we 
move to strong disorder. 

In the 2D case our results 22 on the FM suggest an 
insulating T = state beyond a critical disorder, with 
dp/dT < 0. At the highest disorder, S = 0.08, p(T) in 
the A type AF remains essentially flat down to T = 0, 
while in the G type phase there is still a low temperature 
downturn. While these are finite size results, we argue 
below why there is an intrinsic reason for transport in 
the AF phases to be less sensitive to antisite disorder 
and domain formation. This is related to the nature of 
electronic wavefunctions in these phases. We consider 
the two phases in succession below. 

In the A type phase the core spin order involves diag- 
onal stripes. 'Up-spin' electrons delocalise on down-spin 
stripes which involve one B diagonal and the two adja- 
cent B' lines. 'Down-spin' electrons delocalise on up-spin 
stripes. The down and up stripes share a common B' 
line. The essential feature is that the electronic wave- 
functions are quasi-one dimensional. The introduction of 
antisite disorder leads to two effects: (i) it hinders propa- 
gation along the stripe, and (ii) allows scattering between 
the stripes leading to an 'expansion' of the wavefunction 
in the transverse direction. On its own, the first effect 
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FIG. 7: Colour online: Magnetic correlations with increasing electron density as one traverses the coexistence regime in an 
antisite disordered background. The first panel shows the structural pattern arising from the ASD. Panel 2 shows the spin 
correlations at \i — [Lfm — —1-8, where the ground state is a (domain) ferromagnet. The extreme right is for \i — [iaf — —1-4, 
where the system has only A type AF correlations. As \i increases from \ifm to \iaf the pattern exhibits coexistence of short 
range FM and AF correlations. 



would have suppressed conduction, but the new matrix 
element between stripes allows a transverse pathway for 
delocalisation. In contrast to the FM case there is no 
'confinement' of the wavefunctions to specific domains, 
and the two competing effects above lead to a finite re- 
sistivity (at least a much weaker upturn) even at strong 
disorder. 

In the clean G type phase the system can be viewed 
as two interpenetrating square lattices, one with up spin 
B sites, the other with down spin B sites. The electrons 
delocalise via the B' sites, and each B' site hosts both up 
and down electron states. Electrons in both spin channels 
are delocalised over the whole system in the absence of 
ASD. The presence of antisite domains leads to scattering 
but no confinement of electrons to the domains. Due 
to the inherently 2D character of the G type electronic 
states, in contrast to quasi- ID for A type, the system has 
a lower resistivity. 

We contrast the resistivity of the three phases in Fig. 6 
for the case of maximum disorder that we have studied. 
These are configurations with S ~ 0.08, but, as we have 
noted, with a fairly high degree of local correlation. The 
FM has a clear low T upturn due to the confinement of 
electrons into domains or pathways created by the ASD. 
The A and G type electronic states, in either spin chan- 



0.7 




FIG. 8: Colour online: The variation of electron demsity with 
chemical potential as the FM to A type AF crossover is tra- 
versed in a high ASD sample (S = 0.08). 



nel, are not confined to the magnetic/structural domains, 
and the resistivity remains comparatively lower. The 
high temperature resistivity is determined by spin dis- 
order scattering, depends weakly on carrier density (the 
FM, A and G phases have different n), and is almost 
temperature independent. 

E. Phase coexistence 

In the clean limit, the increase of electron density 
by doping the FM would encounter a window of phase 
separation. A homogeneous state is not allowed for 
0.45 < n < 0.58 and this (idealised) system would break 
up into macroscopic regions having densities n ~ 0.45 
and n ~ 0.58. This pathology is avoided by long range 
Coulomb interactions or quenched disorder. For the anti- 
site disordered configurations that we are considering the 
ASD itself controls the pattern of spatial coexistence. 

Fig. 7 shows the magnetic correlations in a fixed antisite 
background for changing electron density. The leftmost 
panel is the structural pattern, showing the antisite do- 
mains. The next five snapshots correspond to increasing 
chemical potential, /i, and consequently the electron den- 
sity. The first panel in this set is a (domain) ferromag- 
net at n < 0.45, the second shows emergence of stripes 
along with the FM regions. The FM regions shrink and 
the linelike patterns become more prominent in the third 
panel. The fourth and fifth panel complete the evolution, 
with FM correlations completely replaced by stripes (of 
both orientation) as in the bottom of middle column in 
Fig. 3. The evolution of the particle density with /i, and 
the rapid change near the phase separation window, are 
shown in Fig. 8. 

V. CONCLUSIONS 

We have studied the survival of the antiferromagnetic 
double perovskite phases in the presence of spatially cor- 
related antisite disorder. We observe that antisite disor- 
der affects the antiferromagnetic order much less strongly 
than it affects ferromagnetism. For a given structural or- 
der parameter 5, the A type AF structure factor follows 
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D A ~ (1 + 5 2 )/2, in contrast to D F ~ 5 2 in the ferro- 
magnet, while the G type phase follows Dq ~ (l + 6') 2 /4. 
So, despite the possibility of large antisite disorder at the 
high electron doping needed to observe the AF phases, 
there is certainly hope of observing these magnetic struc- 
tures. The AF states are metallic, and the electronic 
wavefunctions in these phases continue to be spatially ex- 
tended even at large disorder. Antisite disorder increases 
the residual resistivity, but, unlike the ferromagnet, we 
did not observe any insulating regime. The field response 



of these AF metals is also fascinating, and will be sepa- 
rately discussed. 
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